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Abstract 

We study in a self-consistent way the impact of the emission of bulk gravi- 
tons on the (homogeneous) cosmology of a three-brane embedded in a five- 
dimensional spacetime. In the low energy regime, we recover the well known 
result that the bulk affects the Friedmann equation only via a radiation-like 
term C/a 4 , called dark or Weyl radiation. By contrast, in the high energy 
regime, we find that the Weyl parameter C is no longer constant but instead 
grows very rapidly as C oc a 4 . As a consequence, the value of C today is not 
a free parameter as usually considered but is a fixed number, which, gener- 
ically, depends only on the number of relativistic degrees of freedom at the 
high/low energy transition. Our estimated amount of Weyl radiation satisfies 
the present nucleosynthesis bounds. 

I. INTRODUCTION 

In the last couple of years, considerable attention has been paid to a new type of models 
with extra dimensions, for which ordinary matter is confined on a three-dimensional space 
called brane. Although a lot of effort has been devoted to investigate the signatures of 
extra-dimensions in gravity tests and collider experiments, these potentially new effects will 
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be inaccessible if the higher- dimensional fundamental mass is too high. Cosmology would 
then be the last hope to detect an extra-dimensional signature, via for instance the Cosmic 
Microwave Background Radiation (CMBR) anisotropics because of the high energies 
reached in the very early universe. 

The simplest viable model of brane cosmology that takes into account the self-gravity of 
the brane is the cosmological extension || of the Randall- Sundrum (RS) setup || and con- 
sists of a ^-symmetric (i.e. mirror symmetric) three-brane, with a tension A and containing 
ordinary cosmological matter with energy density p, embedded in an empty five-dimensional 
spacetime endowed with a (negative) cosmological constant A = —6p 2 . In this context, 
one can solve explicitly |4]], for any equation of state of brane matter, the five- dimensional 
Einstein's equations 



where k 2 = M -3 is the five-dimensional gravitational coupling and where the energy- 
momentum tensor here vanishes everywhere except on the brane. Imposing the RS condition 
H k 2 \ = 6/i, this leads to the generalized Friedmann law 



where H = a/ a is the expansion rate and C is an arbitrary constant. Although the first term 
on the right-hand side, with the identification k 4 A/6 = K 2 p = 8ttG = Mp 2 , agrees with the 
familiar Friedmann equation, the two additional terms induce a deviation from the usual 
law. First, there is a term quadratic in p 0, which dominates the usual linear term in the 
high energy regime, when p ^> A. Second, one finds a radiation-like term proportional to 
the integration constant C, which we will call here the Weyl parameter since it is related to 
the bulk Weyl tensor |J and which is formally analogous to the Schwarzschild mass which 
arises as an integration constant when solving the vacuum Einstein equations with spherical 
symmetry. Indeed, the five-dimensional spacetime metric is Schwarzschild-Anti de Sitter 
(Sch-AdS) which can be written in the manifestly static form || 



Gab + kgAB = k Tab, 



(1) 



18 6b a 4 



(2) 
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ds 2 = -f(r)dT 2 + ^- + r 2 dx 2 , f(r) = „V - ~ 2 , (3) 
j[r) r z 

and the 'radial' coordinate r of the brane can be identified with its scale factor a. When 

the bulk spacetime is empty, C is necessarily constant in time. This can be seen as a 

generalization to five dimensions of Birkoff's theorem M. 

In this setup, one recovers the usual conservation law p + 3H(p + p) — 0. This picture 

is however oversimplified in the sense that, if the assumption of homogeneity is justified on 

average, there exist fluctuations on small scales, which create gravitational waves that can 

escape into the bulk. In the language of particle physics, one would say that interactions 

of brane particles will generate bulk gravitons. If this question has been analysed in detail 

in the context of flat compact extra-dimensions (see e.g. @), it has received so far much 



less attention in the RS context ||10|| , the main difficulty being to take into account the 
back reaction of the graviton flow on the bulk metric. The purpose of the present work is 
precisely to present a treatment that combines self- consistently the emission of gravitons 
with the surrounding spacetime in which our brane-universe is moving. 

II. BRANE IN A RADIATIVE SPACETIME 

In order to model the bulk spacetime metric, we generalize to five dimensions the four- 
dimensional Vaidya's metric [IT|, which is used in general relativity to describe the spacetime 



surrounding a radiating star. Our metric ansatz is given by 

ds 2 = -f(r,v) dv 2 + 2drdv + r 2 dx 2 , f(r, v) = p 2 r 2 - (4) 



If C does not depend on v, the above metric is simply a rewriting of the Sch-AdS metric ([3D 
with the change of coordinate v = T + J dr/f(r). In all cases, v is a null coordinate since the 
spacetime trajectory defined by v — const (and x l = const) is light-like. In Schwarzschild 
spacetime, it would correspond to ingoing radial light rays. Note that, strictly speaking, 
the Vaidya spacetime corresponds to an outgoing radiation flow, whereas in our case, we 
are interested in an ingoing radial flow because the brane, which emits the radiation, is in 



some sense located at the largest radius of spacetime (when getting away from the brane, 
one moves radially 'inwards', as can be deduced from the junction conditions (|8|)). 

In the above spacetime (f|), the trajectory of the brane can be expressed in terms of its 
coordinates v(t) and r(t) as functions of the proper time t (which is also the cosmic time). 
The normalization of the velocity vector (u A = {v, r, 0}, where a dot stands for a derivation 
with respect to t) implies 

r + Vf + r 2 

» = j , (5) 

assuming that the brane is always in the region of spacetime where / > (i.e. outside the 
horizon in the Schwarzschild case). 

Similarly to its four-dimensional counterpart, it is easy to verify that the metric ([§) is a 
solution of Einstein's equations (|l|) with a bulk energy-momentum tensor of the form 

Tab = crk A k B , (6) 

where k A is a (here ingoing) null vector. Our model thus corresponds to a simplified picture, 
where the bulk gravitons are supposed to escape 'radially' from the brane. Choosing to 
normalize it so that Uau a = 1, the only non- vanishing component is k r (= k v ) = 1/v. The 
quantity a represents, for a brane observer, the energy flux of the bulk gravitons and is 
related, via Einstein's equations ([!]), to the variation of the Weyl parameter, according to 
the expression 

dC 2k 2 a 3 /. r~ -\ 2 

-r 3 (r-Jf + r 2 ) . (7) 



dv 3 

The Einstein equations being solved in the bulk, it remains to impose at the brane 
location the Israel junction conditions, which relate the jump of the extrinsic curvature 
tensor for the brane Kab = hf^V ' c n B (n A is the unit vector normal to the brane and 
Hab = 9ab — TiATiB the induced metric) to the brane energy-momentum tensor tab, namely 



[Kab] = k 2 ( tab - \-rh A B ) • (8) 
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Here, tab = (pb + PbjUAUB + PbhAB, where pb = A + p and p& = —A + p are respectively 
the total energy density and pressure in the brane. The ordinary spatial components of the 
junction conditions yields (using ua = {r, —v, 0}) the expression 



-TrPb = , (9) 

o r 

which implies the same unconventional Friedmann equation as in (|2|), with the only difference 
that C is a function of v instead of a constant. The other components (along the time and 
the fifth dimension) yield another relation, which upon using (§) and its time derivative as 
well as (^), can be rewritten as 

Pb + 3-(p b + p b ) = -2ir. (10) 

r 

This equation simply expresses the fact that, from the point of view of the brane, the energy 
loss suffered by the brane exactly corresponds to the bulk radiation flow embodied in (||) 
(note that the right-hand side is very simple here because we have chosen to normalize k A , 
and thus a, with respect to the brane observers). There is a factor 2 because the brane is 
radiating on both sides, into two copies of the same bulk spacetime. 



III. PRODUCTION RATE OF BULK GRAVITONS 

In order to determine quantitatively the energy loss a we now need to evaluate the 
cross-section of the process ip + ip — > G, where ip is a particle confined on the brane (ip 
its antiparticle) and G a bulk graviton, process which appears to be the dominant effect 
in the very early universe. To compute this cross-section, one can ignore the cosmological 
influence since the temperature T of the brane particles is always much bigger than H 
(indeed H ~ T 4 /M 3 in the high energy regime and H ~ T 2 /M p in the low energy regime, 
and we always assume T M) . A Minkowski background in the brane corresponds to the 
five-dimensional Randall- Sundrum solution Starting from the linear perturbations of 
the RS metric, 
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d s 2 = [e- 2 ^ y \ 7]^ + 2Kh^ (x, y)) dx^ dx u + dy 2 , (11) 

in the RS gauge (i.e. h^ a = and d^h au = 0), the gravitons correspond to the decomposition 
into (generalized) Kaluza-Klein (KK) modes, 

V V) = I dmu m (y)<P { ™\x) } (12) 



where the modes u m (y) are given by Wl 

[WJi( m /v) N 2 (me^/n) - N 1 {m/p)J 2 (me^f/j) 

u m (y) = \ -, (13) 

V 2 ^ ^/(J 1 (m//x)) 2 + (iV 1 (m//x)) 2 

and satisfy the normalization / dy e 2 ^ y \ u* m (y)u m i(i)) = 5{m — ml). The h av have been 
defined so that the kinetic part of the graviton action is of the canonical form S^ n — 
— (1/2) / dm J d A x(d(p m ) 2 . The interaction between brane matter (living at y = 0) and the 
bulk gravitons is described by the action 

S int = kJ dmu m (0) J d 4 xr^ 0, (14) 

which means that the effective coupling constant for the 'canonical' KK modes m is 
«|w m (0)|. For the modes m ^> /x, (|T3|) gives us |u m (0)| ~ \ j yfn. 

The interaction Lagrangian (|T4]) gives us the explicit coupling between any bulk KK 
graviton and brane matter, from which one can derive the amplitude for the scattering of 
brane particles leading to a KK emission. At this stage, the calculation is quite analogous 
to the procedure already described in the context of flat extra dimensions [Tj|, the only 
difference being the coupling constant in (|T4|). Using those results, one finds that the spin 
and particle-antiparticle averaged squared amplitude (neglecting the mass of the incoming 
particles) is given here by 

J2\M\ 2 = K 2 \u m (0)\ 2 A S ^, (15) 

o 

with s = (pi + P2) 2 (pi and p 2 being the incoming 4-momenta of the scattering particles), 
and with A s = 2/3 for scalar-scalar scattering, Af — 1 for fermion-fermion scattering, and 
A v = 4 for photon-photon scattering. 
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Going back to cosmology, the brane matter energy density deficit corresponds to the 
energy density created in the form of bulk gravitons, which can be expressed, after integrating 
the Boltzmann equation over momentum space, as 

! + 3tf(, + p) = -/*"/ f^C m , < 16) 

where the collision term associated with the scattering process described above is given by 

d 3 pi d 3 p 2 
(2tt) 3 2Ei (2tt) 3 2E 2 

with the Fermi/Bose distribution functions = l/(e Ei ^ T ± 1), while p m is the four- 
momentum of the created bulk graviton. We will assume that the temperature satisfies 
T ^> /i (this is an excellent approximation in the early universe for a large range of values 
for fi since the only constraint, from gravity experiments, is fi > 10~ 3 eV ). The typical 
energy of gravitons created by brane particles will be of the order of T, so that most of the 
contribution to the energy loss will come from 'heavy' gravitons satisfying m ^> fi. There- 
fore, we simplify the calculation by directly replacing u m (0) by l/y/n in the right-hand side 
of fllED, thus introducing only a slight error due to the very light KK modes (m < fi). One 
can then perform explicitly the integration to find 

"315C(9/2)C(7/2)" 



p + AHp 



g{T)^T\ (18) 



5127T 3 

where we have introduced an effective number of degrees of freedom (relevant for bulk 
graviton production) defined by 

g(T) = ((2/3)0. + Ag v + (l - 2~ 7 / 2 ) (l - 2^) g f ) , (19) 

where g s , gj and g v are the number of relativistic degrees of freedom, at temperature T, for 
particles of respectively spin 0, 1/2 and 1. 



IV. COUPLED COSMOLOGICAL EVOLUTION 



We now combine all the above results in order to get a dynamical system describing 
the coupled evolution of the brane energy density, of the brane scale factor and of the Weyl 
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parameter. 

In the radiation era, on which we focus, the energy density loss rate a can be rewritten 

as 



, = « kV , a= ^ c(9/2)c(7/2) |, (20) 



after use of the relation p = {jt 2 / 30) g*T 4 , where g* = g s + g v + (7/8) <?/ is the total effective 
number of relativistic degrees of freedom in the energy density If all degrees of freedom 
of the standard model are relativistic, = 106.75 and g ~ 166.21 so that a ~ 0.019. 
Introducing the dimensionless quantities p = p/ A, H = H/p, t = pt, the equation governing 
the evolution of the energy density now reads 

§ + AHp=-ap 2 , (21) 
at 

whereas the Friedmann equation takes the form 

H 2 = 2p + p 2 + ^, (22) 

(X 

with the dimensionless Weyl parameter C = C/p 2 . The final equation we need is that 
governing the evolution of the Weyl parameter, which is given by 

^ = 2aa A p 2 {l + p-H). (23) 

It is obtained by writing C = v(dC/dv) and by combining ([?]), (|5]) and the junction condition 
(0) . Note that one must always have H < 1 + p as a consequence of our previous assumption 
/>0. 

The coupled system consisting of (^1|), ( p2|) and ( p3| ) can be analysed both numerically 
and, in some specific regimes, analytically. In the high energy regime p ^> 1, the system 
approaches an asymptotic behaviour characterized by a rapid growth of the Weyl parameter 
due to an abundant production of bulk gravitons: 

C~^—a\ a-t 1 /^. (24) 
4 + a v ' 



One thus sees that the evolution of the scale factor is different from the situation where the 
graviton outflow is ignored, in which case one finds, in the high energy regime, a oc t 1 / 4 . 
However, since a is very small for realistic cases, the power law is essentially the same. 

In the low energy (p -C 1) radiation era, one recovers the conventional evolution a ~ t 1//2 
and the Weyl parameter approaches a constant value. Not surprisingly, this means that, 
during the late conventional cosmological evolution, the production of bulk gravitons can 
be safely ignored. 

However, what is much more remarkable is that the asymptotic value for C is determined 
quantitatively from the evolution of the brane-bulk system in the radiation phase. This 
asymptotic value for C can be estimated analytically. If the present description is valid deep 
enough in the high energy regime, one finds 

PWeyl C a 

p rad 2a 4 p 4 

after the high/low energy transition. It must be stressed that ew depends only on the number 
of relativistic degrees of freedom at the high/low energy transition, and, in particular, is 
independent of the five- dimensional Planck mass. 

The ratio ew, extrapolated at later times (i.e. taking into account the variation of g*), 
is constrained by the number of additional relativistic degrees of freedom allowed during 
nucleosynthesis [|14]], which is usually expressed as the number of additional light neutrino 
species AN U . A typical bound AN U < 1 implies ew ~ 8 x 10~ 2 at nucleosynthesis, whereas 
our model yields the value 

^ 2 x 10 -3 (im.) 1,3 713 (26) 

1, 166.2 1) \ 10.75 ) \g,(T t )) k ' 

where T t is the temperature at the transition between the high and low energy regimes. 

It is rather striking that the Weyl contribution which we have estimated turns out to be 
just below the present observational bounds. Given that forthcoming observations like the 
precise measurement of CMBR anisotropies might improve significantly the above bound 



|15|, this opens the fascinating perspective that the simplest category of brane cosmology 



models could be confronted soon to observational tests at the level of homogeneous cosmol- 
ogy, and not only from cosmological perturbations. 
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